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Abstract. We analyze final-time dependent discrete-time quantum walks in one dimension. We compute asymptotics 
of the return probability of the quantum walk by a path counting approach. Moreover, we discuss a relation between 
04 , the quantum walk and the corresponding final-time dependent classical random walk. 
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1 Introduction 

The discrete-time quantum walk (QW) was first intensively studied by Ambainis et al. [JJ and has been 



widely studied for a decade as a quantum counterpart of the random walk. There are reviews, for example 
Kempe [5] , Kendon [5] , Venegas-Andraca [UJ , Konno [TU] . In this paper we consider the final-time dependent 
discrete-time QW determined by a sequence of 2 x 2 unitary matrices. Particularly, we focus on the return 
probability which is the probability that the walker returns to the origin. As for the work on the return 
probability on QWs, see Stefahak et al. [TS] UHl , for example. 



Mackay et al. [TS] and Ribeiro et al. [TB] studied random unitary sequences and suggested numerically 
that the probability distribution of the QW converges to a binomial distribution by averaging over many 
trials. Their result was proved by Konno [S]. Ribeiro et al. [TB] also reported that periodic sequence is ballistic 
and Fibonacci one is sub-ballistic by numerical evaluations. Machida and Konno [14] presented weak limit 
theorems for time-dependent QWs by the Fourier analysis. Romanelli [17) introduced a time-dependent QW 
and found an analytical expression for the asymptotics of the standard deviation by a continuum limit. This 
expression yields that the parameter a (defined by Eq. (JJ in this paper) determines the degree of the spread 
of the QW for large time. The result is consistent with the numerical result. 
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Related to the return probability, localization is one of important topics in the study of QW. We say 
that a model exhibits localization when a long-time limit of the return probability remains positive. For 
the time-independent and space-homogeneous case, Inui et al. [I] first showed that a three-state QW on 
Z = {0, ±1,±2, . . .} exhibits localization. For the time-independent but space-inhomogeneous case, indeed 
the Hadamard walk with one defect, Konno [T3] proved that localization occurs for a two-state QW on 
Z. Furthermore, Cantero et al. [2 provided explicit expressions for the asymptotic return probabilities for 
a more general setting. In the present paper, we analyze the return probability of the QW, which has 
a connection with Romanelli [T7], and the related QW. Furthermore, we compute the probability of the 
corresponding final-time dependent classical random walk (CRW). 

The rest of the paper is organized as follows. Section 2 treats the definition of the model. In Sect. 3, we 
present main results (Theorems 13.11 and I3.3[) of this paper. Section 4 is devoted to proofs of both theorems. 
The results of the corresponding classical model are given in Sect. 5. We give the proofs of these results in 
Sect. 6 and finally in Sect. 7, we discuss our results. 



2 Definition of the walk 

In this section, we give the definition of a two-state QW with final-time dependence on Z considered here, 
where Z is the set of integers. Let Z + = {1, 2, . . .} be the set of positive integers. For the general setting, 
the time evolution of the walk is determined by a sequence of 2 x 2 unitary matrices, {U n : n € Z + }, where 



U„ = 



(In 



with a n ,b n , c n , d n G C and C is the set of complex numbers. The subscript n means the final time. 

The discrete-time QW is a quantum version of the classical random walk with additional degree of 
freedom called chirality. The chirality takes values left and right, and it means the direction of the motion 
of the walker. At each time step, if the walker has the left chirality, it moves one step to the left, and if it 
has the right chirality, it moves one step to the right. Let define 



\L) = 



\R) = 



where L and R refer to the left and right chirality state, respectively. To define the dynamics of our model, 
we divide U n into two matrices: 



Pn = 







Qn — 



o 



o 



with U n = P n + Q n . The important point is that P n (resp. Q n ) represents that the walker moves to the left 
(resp. right) at position x at each time step. 

For fixed the final time n, we let S^™^/, m) denote the sum of all paths starting from the origin in the 
trajectory consisting of I steps left and to steps right at time k(< n) with I + m = k. For example, when 
n > 3, we have 



^2 "0 Qn-Rn PriQni ^3 (2, 1) Qn-P n ~\~ RnQnRn R n Qn- 



Remark that 2^ (/, m) (1 < k < n) is determined by not Ux, U2, ■ ■ ■ , Uu but only U„. Let be the position 
of our quantum walker at time k. The probability that the walker is in position x at time fc(< n) starting 
from the origin with the initial qubit state ip* is defined by 



P(X k 



where k — I + to, x = —I + m and ip* is a linear combination of \L) and \R) with ||<y9*|| 2 = 1. In particular, 
we take (/?* = T [l/-\/2, i/V%], where T is the transposed operator. Then the probability distribution of the 
walk is symmetric. The following is the important quantity of this paper. 



Pn {0) = P(X n = 0). 
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This is the return probability at the time n. Therefore, £>2n+i(0) = for n > 0. 

In the present paper, we concentrate on the unitary matrix with fixed parameters a,r > 0, which 
was recently introduced and studied by Romanelli |17j : 



U H = 



1-f? 



2n 



(1) 



However, the evolution of his model is different from that of our QW. Because for fixed final time n, his QW 
at time k & {1, 2, . . . , n} is determined by Uj^. On the other hand, that of our QW at time k £ {1,2, . . . ,n} 
is given by which depends only on the final time n. Moreover we consider the following with a, r > 0: 



U K = 



1- ? 



(5)' 



2c> 



The two models defined by the unitary matrices £7^ and U~f seem to be similar, but as we will show in the 
next section, asymptotic behaviors of the return probability are different. Furthermore, if U„ with a = 0, 
then P2n(0) = for any n > 1. On the other hand, if t/j^ with a = 0, then P2n(0) = 1 for any n > 0. For 
a time-independent QW, the Hadamard walk has been extensively investigated by many researchers, which 
is defined by U± (or equivalently, U^) with a = 1 and r = l/v2- In the symmetric Hadamard walk, it is 
known that 

lim np 2n {0) = -, 

n— >oo 7T 

see Konno [121 113] , for example. For the corresponding symmetric classical random walk, 

lim Vnp 2n (0) = —f=. 

n->oo V7T 

3 Our results 

In this section, we give results on asymptotics of the return probability of final-time dependent ID QWs 
determined by and U*£ , respectively. 

Theorem 3.1 For final-time dependent ID QW determined by with a,T > 0, we have 

P2„(0) ~ J (T a (2n) 1 - a ) 2 + J 1 (T a (2n) 1 ~ a ) 2 , 

where P2n(0) is the return probability at the origin at the final time 2n, Jk(x) is the Bessel function of the 
first kind of order k, and f(n) ~ g{n) means f(n)/g(n) — > 1 (n — > oo). 

The proof of this theorem will appear in the next section. By using the result, we have 
Corollary 3.2 (i)lf0<a< 1, then 



lim n 1 - a p 2 „{0) = 



1 



(ii) If a = 1, then 



lim P2n (0) = J (r) 2 + Ji(r) 2 



(Hi) If a > 1, then 



lim n 2 ^ {l-p 2n (0)}= (J) 



2n 
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In fact, this corollary (a ^ 1) corresponds to Eq. (21) in [17] despite the difference of two models. However, 
it is not clear whether the corollary (a — 1) corresponds to his result. Because our result does not have Inn 
term. 

From now on we will give a proof of this corollary. As for part (i), we note that as x — > oo, 

FT HT 

Jq(x) ~ \ — cosfa; — 7r/4), J\(x) ~ \ — sin(x — 7r/4), 
V ttx \ ttx 

see Watson [22]. Thus Jo (a;) 2 + J\(x) 2 ~ 2 /ttx. Then we have the desired conclusion. Part (ii) is a direct 
consequence of the theorem. Concerning part (iii), we see that for small x, 

X X X 

JoO)~l-— , Ji(x) ~ - - — , 

see [22]. So we have Jq{x) 2 + J\{x) 2 ~ 1 — x 2 /4. From this, the proof of part (iii) is complete. 
Next we consider the QW given by U% . 

Theorem 3.3 For time- dependent ID QW determined by U„ with a,r > 0, we have 

2a 



P2 



l(0) ~ (k) i J o(r a (2nY- a r + J 1 (r«(2n) 1 -") 2 } 



The proof of this theorem can be found in the next section. From this theorem, we obtain the following 
result as in Corollary 13.21 So we will omit the proof. 



Corollary 3.4 (i)IfO<a< I, then 



lim n a+1 P2n (0) = - (~ 

n— too 7T \ Z 



(ii) If a = 1, then 



(iii) If a > 1, then 



lim n 2 P2n (0) = (jV {J (r) 2 + Jx(r) 2 }. 



lim n 2 



Corollaries l3.2l and 13.41 show a difference of the asymptotic behavior of the return probability of the two 
models defined by and U„ ■ Indeed, localization occurs for a > 1 case of the model defined by C/jf , but 
there is no localization for any a > in the model defined by J7„ . Furthermore, the rates of convergence of 
the return probability of the two models are obtained. 



4 Proofs of Theorems 13.11 and 13.3 



First we will prove Theorem 13.11 Let x„ = (n/r) 20 and A n — \b n \/\a n \. So \a n \ — l/^/x^ and A n 
\Jx n — \ ~ \fx~^. We assume that n is even. We begin with 

"/2 , ,„ v 2 n/2-l 2 

Kn = \ an r 2 x: (-Aiy- 1 r 2 _- 1 ) = Kr 2 £ a^r ( n/2 ~ 1 



. T — 1 / ^— ' V 7 

7=1 x ' 7 7=0 v ' 

In order to compute the asymptotic behavior of K n as n —¥ oo, we introduce 
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Let [z n ] (f(z)) denote the coefficient of z n in the formal power series f(z) = Y^=o fnZ n , so that [z n ] (f(z)) = 
f n . Then we should remark that 

K n =[z°]($ n (z)). (2) 
So we will calculate the asymptotic behavior of & n (z) as n — > oo. Put w = z + z^ 1 . Then we have 

*„(*) = { Kl 2 (1 + iwA n - A*) Y' 2 - 1 ~ (-1)"/ 2 - 1 (1 - lw /^-) n/2 - 1 . 
Noting that (ra/2 - l)/y/x^ ~ T a n 1 ^ a /2, we obtain 

~ cxp (-ir^ 1 - (z + z- 1 ) /2) . (3) 

On the other hand, the generating function of the modified Bessel function is given by 

oo 

exp(u(z + z- 1 )/2)= I m{u)z m , (4) 



771— — OO 



see By Eqs. ©, and Q, we have 

[z°] (*„(*)) ~ {-l) n / 2 - 1 I Q {-iT a n 1 - a ). (5) 

From iz) — Jo{z), we get 

[z°] ($ n (z)) ~ (-l)»/ 2 - 1 Jo(r a n 1 - c ). 

Therefore Eq. © yields 

~ (-l) n / 2 - 1 Jo(r a n 1 - a ). (6) 

Next we consider 

In order to compute the asymptotic behavior of L n as n — > oo, we introduce 

$„(z) = |a„|"- 2 (1 + iA n z) n/2 (1 + iAnZ- 1 ) 71 ' 2 ' 1 . 
Then we should remark that 

£ * = -^{w( $ »w)}- (7) 

Noting that $„(z) = (1 + iA n z) $> n (z), we see that 

[z] (*„(*)) = [z] ($„(*)) + iA n [z°) ($„(*)) . (8) 

By Eqs. ©, and (gj, we have 

[z] ($„(z))^(-l)"/ 2 - 1 / 1 (- l r Q n 1 - Q ). 
From Ii(-iz) = (-i)Ji(z) (see 22), we get 

[z] ($„(*)) ~ (-lr/^HyxCrV-). (9) 
By using Eqs. ©, ©, Q, ©, ©, and A n ~ we obtain 

Here we use the following result given by Konno [3 [5] . 
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Lemma 4.1 When n is even, we have 



P(X n = 0) = \a n \ 2 ^J2J2(- A 

7=1 8=1 



2 s 7 +i5 / n/2 — 1 \ (n/2 — 1 



7-1 / V 5- 1 



2 1 n (1 l\ 1 



2/ 7 <5 2 V7 <V 1 - l a 



2 



kr(-<r u o ) i« - + ( 1 _ K|2 ) 



1 > ~ 2 



n f 5 



w/iere A„ = |6„|/|a n |. 

By using this lemma, noting that n is even, we see that 



P(X n = 0) ~ Xn [{^) 2 Ll - nL nKn + ( T -i 7 -) 



K 



■>;, I - 1 )" t{Jv- -=) 2 - 2 ( J - -L) •/„ ■ fl ■ —)./,; 



~ f 2 4- J 2 

where J fc = ifT-n 1 -") for k = 0, 1. So 

P(X 2n = 0) ~ J (r«(2n) 1 - Q ) 2 + J 1 (r«(2n) 1 - Q ) 2 

completes the proof of Theorem 13.11 

From now on we will prove Theorem 13.31 In this case, we should note that \a n \ = y/l — l/x n and 
A n = \b n \/\a n \ ~ 1/y/x^. Similarly, we see that 

$„(z) = | a „|"- 2 (1 + iA n z) n/2 - x (1 + lAnZ- 1 ) 11 ' 2 - 1 
= {\a n \ 2 (l + lW A n -Al)} n/2 - 1 



Thus we have 



Then 



(1 + lW/y/x n ) 



$ n (z) ~ exp (j^n 1 " (z + z" 1 ) /2) 



n 

By using Lemma l4.il noting that n is even, we see that 

P(X n = 0) ~ ( 1 - ^ | (|) 2 L 2 - ni„X„ + x„K 2 



~ -J {(Jo + V^n^l) 2 - 2 (J + \/x^Jl) Jo + Jo } 

where Jk = Jfe(r Q n 1_Q ) for fc = 0, 1. Therefore the desired conclusion is obtained. 
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5 Final-time dependent classical random walk 

In this section we study the classical counterpart of our final-time dependent QW. We first consider the 
evolution of the time- independent (i.e., usual) CRW on the line is given by 

P( the particle moves one unit to the left) 

p, if the previous step was to the left, 



if the previous step was to the right, 



and 



P( the particle moves one unit to the right) 

1 — p, if the previous step was to the left, 

q. if the previous step was to the right. 

This model is also called the correlated random walk (see [TT], for example). When p — q, the probability 
that the particle moves one unit in the same direction as the last step is p, and the probability that the 
particle moves one unit in the opposite direction as the last step is 1 — p. In particular, if p = q = 1/2, then 
the walk is equivalent to the well-known symmetric (non-correlated) random walk, i.e., the particle moves 
at each step either one unit to the left with probability 1/2 or one unit to the right with probability 1/2. 
The directions of different steps are independent of each other. 

Here we present the definition of the final-time dependent CRW. For the general final-time dependent 
setting, the time evolution of the CRW is determined by a sequence of 2x 2 transition matrices, {V n : n g Z + }, 
where 



d„ 



with a n ,b n , c„ , d n e [0, 1] and a n + Cn — b n + d n — 1. The subscript n indicates the final time. In particular, 
we investigate the following model corresponding to the QW determined by with (3, 9 > 0: 



V R = 



. n. J 

-e\0 



1 



Moreover we consider the following model corresponding to the QW determined by with /3, 9 > 0: 

1 



V 



K 



1 



In this paper, we choose ip* = T [l/2, 1/2] as the initial distribution of the final-time dependent CRWs. We 
first give a result on the return probability of the CRW defined by V^. 

Theorem 5.1 For final-time dependent ID CRW determined by with (3, 9 > 0, we have 

p 2n (0) ~ cxp (-^(2n) 1 -^) {I {9 f \2n) 1 -^) + 7 1 (^(2n) 1 -' 3 )} , 

where P2n{0) is the return probability at the origin at the final time 2n and Ik(%) is the modified Bessel 
function of the first kind of order k. 

The proof will appear in the next section. By this theorem, we have 
Corollary 5.2 (i) If < (3 < 1, then 

'2 



lim n 

n— >oc 



<'-«^„ ( o ) = -L(! 



0/2 
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(ii) If (3 = 1, then 

lim p 2n {0) = e- e {I Q (6) + h{0)} 

(iii) If p > 1, then 



lim n^ 1 {l-P2n(0)}= f|Y 

n— >oo \ z / 



From now on a proof of the corollary will be given. Concerning part (i), it is sufficient to note that as x — > oo, 

e x / 1 \ e x ( 3 
/ (x) == 1 + — , 7 x (a;) == 1 

see [35]. Part (ii) is trivial. As for part (iii), we observe that for small x, 

X 3j X 

I (x)~l + — , J 1 ( a; )~- + — ) 

see [22]. So we get Io(a;) + Ji(2;) ~ l + x/2 and then e~ K {/o(^) + h(x)} ~ l — x/2. Thus we have the desired 
conclusion. 

Next we present a result on the return probability of the CRW given by Vjj* . 
Theorem 5.3 For final-time dependent ID CRW determined by with /3,t > 0, we have 

/ ft \ 13 

P2n(0) ~ (—) exp (-^(2n) 1 ^) {/ (^(2n) 1 -' 3 ) + h{6^ {2n) 1 ^)} . 

The proof can be seen in the next section. From this theorem, we obtain the following result as in Corollary 
15.21 So we will omit the proof. 

Corollary 5.4 (i)If0<f3< 1, then 



(ii) If ft = 1 , then 



(iii) If (3 > 1, then 



Urn nV+Wp 2n (0) = 4= f? 

n— >oo -i/7T \ Z 



lim np 2 „(0) = ^— {I o (0) + h(0)}. 

71— >00 Z 



lim n 13 p 2n (0) 



6 Proofs of Theorems 15.11 and 15.3 



First we will give a proof of Theorem 15.11 Let y n — (n/9)° and B n — b n /a n . Note that a n = l/y n and 
B n = y n — 1. We assume that n is even. We begin with 

"/ 2 / /o 1\ 2 n/2-1 / . N 2 

*<r> = «r 2 E { ! r _\ ) = <- 2 E B » 7 ( n/ 7 " 

In order to compute the asymptotic behavior of as n —¥ oo, we introduce 



* n (z) = ar 2 (1 + B^f 12 - 1 (1 + BnZ^ 1 ) 



_ 1 \n/2-l 



Note that 

K^ = [z°] (*„(*)). 

So we will calculate the asymptotic behavior of <$>n( z ) as n — >■ oo. Then we have 

*„(*) = {a 2 n (1 + wB n + B 2 n ) Y' 2 - 1 ~ {1 + (w - 2)/y n } n/2 - 1 . 
Noting that (n/2 - l)/y n - 6» /3 ti 1 " /3 /2, we obtain 

OO 

~ exp (-^n 1 -' 3 ) exp (^n 1 ^ (z + z" 1 ) /2) - exp (-^n 1 ^) ]T J^n 1 "' 3 ); 

m— oo 

Then wc have 

[z°] ($„(*)) ~ exp (-^n 1 "") I (^n 1 -' 3 ). 

Therefore 

ffW~eKp(-0V-*)lo(0V-''). 



Next we consider 



/ 2 1 /_ /n ,\ 2 _„_ 2 n/2-1 

n 

7 v 7-1 y ^ " n V7 + iy v 7 



^l/n/2-lV 2ar 2 U d27+1 / n/2 \ (n/2 - 1 



L (c) n-2 

7—1 ' ' (t 7—0 



^ 71 7 1-1 " nB„ ^ n 



In order to compute the asymptotic behavior of L$ as n — > oo, we introduce 

$„(z) = «r 2 (i + B nZ ) n/2 (i + B^y' 2 - 1 . 

Then we should remark that 

Noting that $,i(z) = (1 + B n z) &n(z), we see that 

[z] ($„(*)) - [z] ($„(«)) + B n [z°] ($„(z)) . 

On the other hand, from Eq. pU|) we have 

[z] ($„(«)) - exp (-0V^) h^n 1 -?). 
Therefore, by Eqs. dTTJ), ([11]), ([13]), and ([11]), we obtain 



~ | exp (-0V^) {/ O (0V^) + ^J^} 

„ 2 exp ( _^ n i-^ ( /o( 0V-/*) + M^!) 



Here we use the following fact given by 
Lemma 6.1 When n is even, we have 

P(Y n = 0) = < /2 ~ 2 < /2 ~ 2 ||(a„ + d n )b n c n L^ + X -{a n c n + b n d n )(a n d n - b n c n )K^ 



By using this lemma, noting that n is even, we see that 

P(Y n = 0) ~ y n { f 1 - - + 4) ?4 C) + f-1 + 1 4) 

IV Vn VnJ 2 V Vn Vl) 

~ exp {/o^n 1 "^) + /^n 1 ^)} . 

Hence, the proof of Theorem 15. 11 is complete. 

Now we will move to prove Theorem 15.31 In this case, we should note that a„ = 1 — l/y n and B n = 
K/a n = l/(y„ - 1) — l/y n . Similarly, we see that 



Thus we have 



Then 



$„(z) = <- 2 (1 + B n zf l2 - X (1 + Bn*- 1 )" 72 " 1 
= {al(l + wB n + Bl)} n/2 - 1 
~ (1 + (w - 2)/y n ) n/2 - 1 . 

$„(«) - exp (-O^n 1 -?) exp (flV"' 3 (z + z' 1 ) /2) 



~ exp /o^n 1 ^), 4 C ) ~ exp (-^n 1 "") - {lo^n 1 -") + ynh^n 1 ^)} 

By using Lemma 16.11 again, we see that 

P(Yn = 0) ~ f 1 - Pj \ {^i C) + (Vn ~ 2)ifW} 

exp (-^n 1 "") {/o + y„/i + (y„ - 2)J } 



exp 



where Ife = I^O^n 1 " 13 ) for fc = 0, 1. Therefore we have the desired conclusion. 



7 Discussion 

In the present paper, we have the asymptotic behavior of the return probability of the final-time dependent 
QWs on the line defined by and tTjf and the corresponding final-time dependent CRWs on the line 
defined by Vj^ and ■ Under the condition a = j3, the rate of convergence of return probability of the QW 
is square of that of the CRW. The QWs (resp. CRWs) defined by with a > 1 (resp. V r f with /3 > 1) 
exhibit localization. On the other hand, the QWs (resp. CRWs) defined by Ui^ with any a (resp. with 
any /3) do not exhibit localization. For the CRW defined by Vjf with (3 = 1, the corresponding weak limit 
theorem was obtained by [II] . Furthermore, the asymptotics of the QWs (resp. CRWs) given by with 
< a < 1 (resp. < (3 < 1) correspond to weak limit theorems, (i.e. Corollary 3 (resp. (A. 7))) appeared in 
Chisaki et al. [3]. 

Acknowledgment. The second author thanks F. Alberto Griinbaum for useful comments on an early 
version of this paper. This work was partially supported by the Grant-in-Aid for Scientific Research (C) of 
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